Electron spin relaxation in p-type GaAs quantum wells 
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We investigate electron spin relaxation in p-type GaAs quantum wells from a fully microscopic 
kinetic spin Bloch equation approach, with all the relevant scatterings, such as the electron-impurity, 
electron-phonon, electron-electron Coulomb, electron-hole Coulomb and electron-hole exchange (the 
Bir-Aronov-Pikus mechanism) scatterings explicitly included. Via this approach, we examine the 
relative importance of the D'yakonov-Perel' and Bir-Aronov-Pikus mechanisms in wide ranges of 
temperature, hole density, excitation density and impurity density, and present a phase-diagram-like 
picture showing the parameter regime where the D'yakonov-Perel' or Bir-Aronov-Pikus mechanism 
is more important. It is discovered that in the impurity-free case the temperature regime where 
the Bir-Aronov-Pikus mechanism is more efficient than the D'yakonov-Perel' one is around the hole 
Fermi temperature for high hole density, regardless of excitation density. However, in the high 
impurity density case with the impurity density being identical to the hole density, this regime is 
roughly from the electron Fermi temperature to the hole Fermi temperature. Moreover, we predict 
that for the impurity-free case, in the regime where the D'yakonov-Perel' mechanism dominates 
the spin relaxation at all temperatures, the temperature dependence of the spin relaxation time 
presents a peak around the hole Fermi temperature, which originates from the electron-hole Coulomb 
scattering. We also predict that at low temperature, the hole-density dependence of the electron 
spin relaxation time exhibits a double-peak structure in the impurity-free case, whereas first a peak 
and then a vaUey in the case of identical impurity and hole densities. These intriguing behaviors 
are due to the contribution from holes in high subbands. 

PACS numbers: 72.25.Rb, 67.30.hj, 71.10.-w 



I. INTRODUCTION 

In recent years, much attention has been devoted to 
semiconductor spintronics both theoretically and exper- 
imentally due to the potential application of spin-based 
devicesjii^ In order to manipulate the spin relaxation 
such that the information is well preserved before re- 
quired operations are completed, it is crucial to gain 
a thorough understanding of spin relaxations. In p- 
doped III-V semiconductors, the main electron spin re- 
laxation mechanisms have been recognized as>i the Bir- 
Aronov-Pikus (BAP) mechanism^ and the D'yakonov- 
Perel' (DP) mechanism.^ In the DP mechanism, elec- 
tron spins decay due to their precessions around the 
momentum-dependent spin-orbit fields (inhomogeneous 
broadening)^ during the free flight between adjacent scat- 
tering events. In the BAP mechanism, spin relaxes due 
to spin- flip caused by exchange interaction with holes. 
It was believed that in p-doped bulk samples the BAP 
mechanism dominates the spin relaxation process at high 
doping density and low temperature, whereas the DP 
mechanism is more important at low doping density 
and high temperaturej^i^ i^'^'^" In two-dimensional sys- 
tem, Maialleii calculated the spin relaxation time (SRT) 
due to these two mechanisms at zero temperature by us- 
ing the single-particle approach and showed that these 
two SRTs have nearly the same order of magnitude. How- 
ever, as pointed out by Zhou and Wu lately;^ there 
are some common problems in the previous literature: 
The SRT due to the BAP mechanism was calculated 



based on the elastic scattering approximation, which is 
invalid at low temperature due to the omission of the 
Pauli blocking. Also, the investigation of the SRT due 
to the DP mechanism was also inadequate because the 
Coulomb scattering is not included in the frame of the 
single-particle theory. 

Zhou and Wu applied the fully microscopic ki- 
netic spin Bloch equation (KSBE) approach^-^"^ to 
investigate the spin relaxation in p-type GaAs quan- 
tum wellsji^ The KSBE approach has achieved good 
success in the study of the spin dynamics in semi- 
conductors, where not only the results are in good 
agreement with the previous experiments, but also 
many predictions have been confirmed by the latest 
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Via this approach, they explicitly included all the rel- 
evant scatterings and obtained the accurate SRT due 
to these two mechanisms. It was found that the 
BAP mechanism is always less efficient than the DP 
mechanism for moderate and high excitation densities 
where A'ex ^ O.liV/i [Ncx (Nh) is the excitation (hole) 
density] , in contrast to the common belief in the previous 
literatureiiiiiSiiSiiiS This claim has very recently been 
confirmed experimentally by Yang et al.^ Moreover, 
a similar conclusion was also obtained in bulk GaAs 
later 

However, for very low excitation density where the 
Pauli blocking of electrons is negligible, for high hole den- 
sity where the contribution from the high subbands or 
different hole bands becomes significant and/or for high 
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impurity density where the spin relaxation due to the DP 
mechanism is suppressed, whether the BAP mechanism 
can be more efficient is still questionable. In the present 
work, we extend the KBSEs to include both the lowest 
subband of light-hole (LH) and the lowest two subbands 
of heavy- hole (HH) , and compare the relative importance 
of the DP and BAP mechanisms in wider ranges of tem- 
perature, hole density, excitation density and impurity 
density. We present a "phase-diagram-like" picture indi- 
cating the dominant spin relaxation mechanism. In the 
case with no impurity and high excitation density, our 
results show that the BAP mechanism is unimportant 
at low temperature, in consistence with Ref. Il2l . Nev- 
ertheless, since more hole subbands and bands are in- 
cluded in our model, we are able to discuss the case with 
higher hole density. We find that the BAP mechanism 
can surpass the DP mechanism at high temperature for 
sufficiently high hole density. In the case with no im- 
purity and low excitation density, the BAP mechanism 
can surpass the DP mechanism for wider hole-density 
and temperature ranges. Moreover, we also find that in 
both cases above, the regime where the BAP mechanism 
is more efficient is always around the hole Fermi tem- 
perature for high hole density, regardless of excitation 
density. However, in the high impurity density case with 
the impurity density being identical to the hole density, 
the behavior is very different from the impurity-free case: 
the regime of hole density where the BAP mechanism is 
more efficient becomes larger, and the regime of temper- 
ature becomes wider, ranging roughly from the electron 
Fermi temperature to the hole Fermi temperature. We 
also show that the multi-hole-subband effect leads to a 
very intriguing hole-density dependence of SRT at low 
temperature. 

This paper is organized as follows. In Sec. II, we set 
up the KSBEs. In Sec. Ill, we compare the relative im- 
portance of the BAP and DP mechanisms in different pa- 
rameter regimes and investigate the multi-hole-subband 
effect. We conclude in Sec. IV. A comparison of the cal- 
culation from the KSBEs with the experimental data of 
a p-type GaAs quantum well is given in the Appendix. 



II. KINETIC SPIN BLOCK EQUATIONS 

We investigate a p-type (001) GaAs quantum well of 
width a with its growth direction along the z-axis. The 
width is assumed to be small enough so that only the low- 
est subband of electron, the lowest two subbands of HH 
and the lowest subband of LH are relevant for the electron 
and hole densities we discuss. The envelope functions 
of the relevant subbands are calculated under the finite- 
well-depth assumptionJi2.iii The barrier layer is chosen 
to be Alo.4Gao.6As where the barrier heights of electron 
and hole are 328 and 177 meV, respectively.— We focus 
on the metallic regime where most of the carriers are in 
extended states. Since the hole spins relax very rapidly 
(only several picoseconds), we assume that the hole sys- 



tem is always in the equilibrium. 

Via the nonequilibrium Green function method,— we 
construct the KSBEs as foUowsj^^ii^ 



dtPW — dtPk\coh + dtPklscat, 



(1) 



with pk representing the electron single-particle density 
matrix with a two-dimensional momentum k = (fc^,, fcj,), 
whose diagonal and off-diagonal elements describe the 
electron distribution function and spin coherence respec- 
tively. The coherent term can be written as (ft = 1 
throughout this paper) 



dtPkl 



coh 



EHF(k), Pk 



(2) 



in which [A, B] = AB — BA is the commutator. h(k) rep- 
resents the spin-orbit coupling of electrons composed of 
the Dresselhaus^ and Rashbai^ terms. For GaAs quan- 
tum wells, the Dresselhaus term is dominanti^ and 



(3) 



h(k) = 27D - {ki)), ky{{ki) ~ ki), 



where (fc^) stands for the average of the operator 
— {d/dzY over the state of the lowest subband of electron, 
and 7d = 8.6 eV-A^ denotes the Dresselhaus spin-orbit 
coupling coefficient. ^^'^^ SHF(k) is the effective magnetic 
field from the Coulomb Hartree-Fock contribution.^^ For 
the screened Coulomb potential, the screening from elec- 
trons and holes is calculated under the random phase 
approximationji^i^i The scattering term 9t/5k|scat con- 
sists of the electron-impurity, electron-phonon, electron- 
electron Coulomb, electron-hole Coulomb, and electron- 
hole exchange scatterings. The expressions of these scat- 
terings are given in detail in Ref. 12. Here we just extend 
the electron-hole Coulomb and exchange scatterings to 
the multi-hole-subband case. The expression of electron- 
hole Coulomb scattering is still similar to that in Ref. [l^ . 
The complete electron-hole exchange scattering term is 
written as 

9t/5k|BAP = -7rXl'^(^k-q - Ck + Ek'a - ek'-q,A') 

k'qAA' 

T,= ± 

X |r,V(k + k' - q)|2 k,p>_qs_,p< (1 - ft>^x)ft'-^.x' 



^'nPk-q^-vPkft\xi^ /k'-q,A') 



h.C. 



(4) 



Here p^ ~ 1 — pk and p^ ~ pk are the electron 
density matrices; s± = s^: ± isy are the electron spin 
ladder operators. A — HH^"', LH*-"-* with the super- 
script being the subband index of hole, f^y^ is the 
hole distribution on the Ath hole band. The matrix 
comes from the long-range term of the electron- 
hole exchange interaction Hamiltonian and can be writ- 
ten as r± = where A£;lt is the 
longitudinal-transverse splitting in bulk; |03d(O)P = 
l/(7raQ) with ao being the exciton Bohr radius; and 



{— (M )''') are operators in hole spin space. The 
matrix M~ is given bj*^ (in the order of ||}*-^\ 
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(5) 

where K = k + k' — q is the center-of-mass momentum 
of the electron-hole pair with K± = ztiKy. The form 
factors can be written as 



dqz 



with 



/aa'(<Z.) = / dzdz'Uz')Q {z')e 



27r 



(6) 



(7) 

In Eq. ([5]), it is seen that most of the nonzero elements 
in matrix M~ contain K^. or K±, and thus the mag- 
nitudes of these terms increase with increasing K. The 
only exception is M^i i , where the K dependence is only 

from the form factor F^m- Consequently the magnitude 
of -M~i 1 decreases with K. This K dependence con- 

2 ' 2 

tributes to an intriguing hole density dependence of spin 
relaxation to be addressed in this work. 

TABLE I: Material parameters used in the calculation 
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-0.44 




m* 
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"^LH.II 
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Ko 


12.9 






10.8 
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5.31 X 10^ 


kg/m^ 


ei4 


1.41 X 10^ V/m 


Vst 


2.48 X 10^ 


m/s 


Vsl 


5.29 X 10^ m/s 




8.5 eV 






35.4 meV 


Aso 


0.341 eV 




E, 


1.55 eV 


ABlt 


0.08 meV 




do 
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III. RESULTS AND DISCUSSIONS 

By numerically solving the KSBEs with all the scat- 
terings explicitly included, one is able to obtain the SRT 
from the temporal evolution of the electron spin polariza- 
tion along the z-axis. We choose initial spin polarization 
P — 4 %4£ and well width a = 10 nm, external mag- 
netic field B = unless otherwise specified. The other 
material parameters are listed in Table HP^i^'^ 
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FIG. 1: (Color online) Ratio of the SRT due to the BAP 
mechanism to that due to the DP mechanism, tbap/tdp, as 
function of temperature and hole density with (a) Ni — 0, 
Afex = 10" cm-2; (b) iV, = 0, A^cx = 10« cm'^; (c) = Nh, 
iVcx = 10" cm-2; (d) Ni = N^, iVex = 10^ cm'^ The black 
dashed curves indicate the cases satisfying tbap/tdp ~ 1- 
Note the smaller the ratio tbap/tdp is, the more important 
the BAP mechanism becomes. The yellow solid curves indi- 
cate the cases satisfying df^^lN-^^id) + N-^-^(2)]/df_i,^Nh — 0.1. 
In the regime above the yellow curve the multi-hole-subband 
effect becomes significant. 



A. Comparison of the BAP and DP mechanisms 

We first examine the relative importance of the BAP 
and DP mechanisms for different parameters in p-type 
GaAs quantum wells. In Fig. [1] the ratio of the SRT due 
to the BAP mechanism to that due to the DP mechanism 
is plotted as function of temperature and hole density in 
the cases with no/high impurity and low/high excita- 
tion densities. From this figure, one can recognize the 
parameter regime where the DP or BAP mechanism is 
more important. It is also shown that the multi-hole- 
subband effect becomes significant for high temperature 
and/or high hole density (the regime above the yellow 
solid curve). Here and hereafter, the multi-hole-subband 
refers to either the high HH subband or the LH subband. 
Although the multi-hole-subband effect has important ef- 
fect on electron spin relaxation in the relevant regime, the 
main physics is still the same as that in the single-hole- 
subband model. Therefore, in this subsection, we first 
discuss the general behavior about how the relative im- 
portance of the BAP and DP mechanisms is influenced 
by the temperature, hole density, excitation density and 
impurity density, which is analogous in both the multi- 
hole-subband and single-hole-subband models. We then 
investigate the special features from the contribution of 
high hole subbands in next subsection. 

In the case with no impurity and high excitation den- 
sity [Fig. [Ha)], our results are consistent with Ref. [l^: 
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FIG. 2; (Color online) SRTs due to the DP and BAP mech- 
anisms, the total SRT together with the ratio tbap/tdp 
vs. temperature T for A'ox = 10^ cm~^ (curves with •), 
3 X 10^" cm~^ (curves with A) and 10^^ cm~^ (curves with 
■ ) with hole density A''^ = 5 x 10^^ cm~^ and impurity den- 
sities (a) Ni = and (b) Ni = A'^^. The electron Fermi 
temperatures for those excitation densities are T§ — 0.41, 
12.4 and 41.5 K, respectively. The hole Fermi temperature is 
Tp = 124 K. Note the scale of tbap/tdp is on the right-hand 
side of the frame. The multi-hole-subband effect is taken into 
account in the calculation. 



i.e., the BAP mechanism is unimportant at low temper- 
ature, which is in stark contrast with the common behef 
in the hteraturejiiii^iiSiiiS Moreover, since we extend the 
scope of our investigation to higher hole density by in- 
cluding more hole subbands in our model, it is discovered 
that the BAP mechanism can surpass the DP mechanism 
in the regime with high temperature and sufficiently high 
hole density (the regime embraced by the black dashed 
curve). 

In the case with no impurity and low excitation density 
[Fig. [IJb)], one can see that the regime where the BAP 
mechanism surpasses the DP mechanism becomes larger. 
The underlying physics is shown in Fig. [2ja). It is seen 
that the SRTs due to the BAP and DP mechanisms both 
decrease with increasing excitation density {Nex=Ne), 
but the amplitude of the latter is much larger than the 



former. The decrease of tdp comes from the increase of 
the inhomogeneous broadening (|hkp) oc Nexr^^ and 
the decrease of tbap is mainly from the increase of the 
average electron velocity (vk) oc N^^^M' Moreover, the 
increase of the Pauli blocking of electrons can partially 
compensate the effect of the increase of {vk)}"^ Con- 
sequently, Tbap decreases with N^x much more slowly 
than Tdp and the relative importance of tbap is en- 
hanced for lower excitation density. It is also noted that 
when the electron system is in the nondegenerate regime 
(T > Tp = E^/k-Qjj the inhomogeneous broadening and 
{vk) is not sensitive to A^cx- Thus the ratio tbap/tdp 
changes little with the excitation density. 

By comparing Fig. [T{a) and (b), it is seen that the 
regimes where the BAP mechanism is more efficient in 
both cases are around the hole Fermi temperature Tp = 
Ep/ks for high hole density. Here Ep represents the 
Fermi energy of hole at zero temperature calculated with 
the HH(i) , LH(i' and HH^^) subbands included. A typical 
case is shown in Fig. [2ja) for Nh = 5 x 10^^ cm~^4i It is 
shown that the ratio tbap/t'dp ffi'st decreases and then 
increases with increasing The minimum is around 
Tp = 124 K, regardless of excitation density. The under- 
lying physics is as follows. On one hand, the SRT due 
to the DP mechanism first increases and then decreases 
with T and the peak appears around the hole Fermi 
temperature. This is because the electron-hole Coulomb 
scattering, which dominates the momentum scattering, 
increases with increasing temperature in the degenerate 
regime (T < Tp) and decreases with T in the nondegen- 
erate regime (T > Tp), similar to the electron-electron 
Coulomb scatteringJIi^iil On the other hand, the SRT 
due to the BAP mechanism first decreases rapidly and 
then slowly with T. The decrease of tbap is mainly from 
the decrease of the Pauli blocking of holes and the in- 
crease of the matrix elements in Eq. ([5]).^^ In high tem- 
perature (nondegenerate) regime, the Pauli blocking be- 
comes very weak, and thus tbap decreases slowly with 
T. Under the combined effect of these two mechanisms, 
the valley in the ratio tbap/^dp appears around Tp. 

Moreover, we also show that in the regime where the 
DP mechanism is dominant at all temperatures, e.g., the 
high excitation density case [the curves with squares in 
Fig. Ufa)] , the total SRT shows a peak around the hole 
Fermi temperature. This temperature dependence is sim- 
ilar to the peak first predicted theoretically and then con- 
firmed experimentally in n-type samples J^^^iii? The only 
difference is that the peak in the previous work comes 
from the electron-electron Coulomb scattering and hence 
appears around the electron Fermi temperature, whereas 
the peak here originates from the electron-hole Coulomb 
scattering and thus appears around the hole Fermi tem- 
perature. 

Then we turn to the case of high impurity density 
with Ni = Nh [Fig. (TJc) and (d)]. In this case, the 
regime where the BAP mechanism is more important be- 
comes larger than that in the impurity-free case. The 
scenario is that the higher impurity density strengthens 
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the electron-impurity scattering and suppresses the DP 
mechanism, consequently enhances the relative impor- 
tance of the BAP mechanism. Interestingly, it is also 
seen that the temperature regime where the BAP mech- 
anism surpasses the DP mechanism in this case is very 
different from that in the impurity-free case. This regime 
is roughly from the electron Fermi temperature to the 
hole Fermi temperature for high hole density. To explore 
the underlying physics, we plot the SRTs due to the BAP 
and DP mechanisms in Fig.[2{b) for = 5x 10^^ cm~^. 
It is seen that the SRT due to the DP mechanism first de- 
creases slowly and then rapidly with temperature. This 
is because the electron-impurity scattering, which domi- 
nates the momentum scattering, has a very weak temper- 
ature dependence. Thus the temperature dependence of 
tdp is mainly determined by the inhomogeneous broad- 
ening from the spin-orbit coupling. It is also noted that 
the inhomogeneous broadening changes little with tem- 
perature when T < Tp, hence r^p varies with T very 
mildly at low temperature. On the contrary, as men- 
tioned above, the SRT due to the BAP mechanism first 
decreases rapidly and then slowly with temperature. As 
a result, the temperature dependence of tbap/tdp can 
be easily understood. When T < Tp, tdf decreases with 
T slower than tbap, thus the ratio decreases with T. 
In the case with T > Tp, tdf decreases with T faster 
than Tbap, hence the ratio increases with T. The ratio 
Tbap /tdp varies mildly when temperature varies from Tp 
to Tp. Consequently, when hole density is high enough, 
the BAP mechanism can surpass the DP mechanism in 
the temperature regime between these two temperatures. 
In particular, in the case with high impurity and very 
low electron excitation densities [Fig. [T]Jd)], the electron 
Fermi temperature (0.41 K) is much lower than the lowest 
temperature (5 K) of our computation and the hole Fermi 
temperature is close to the highest temperature (300 K) 
of our computation. As a result, the BAP mechanism 
dominates the spin relaxation in the whole temperature 
regime of our investigation. 



We stress that the different behaviors in the impurity- 
free and high impurity density cases originate from the 
different dominant momentum scatterings: the electron- 
hole Coulomb scattering in the impurity-free case and the 
electron-impurity scattering in the high impurity density 
case. The different dominant scatterings lead to the dif- 
ferent temperature dependences of tdp, and hence the 
different behaviors of the ratio tbap/tdp- In the case 
with moderate impurity density, these two scatterings 
both contribute to the DP spin relaxation, thus the trend 
of the temperature dependence of tdp is between those 
in the impurity- free and high impurity density cases. As 
a result, the temperature regime where the BAP mech- 
anism is more efficient than the DP mechanism is from 
some temperature between the electron and hole Fermi 
temperatures to the hole Fermi temperature. 



B. Multi-hole-subband effect 

Now we investigate the multi-hole-subband effect on 
the spin relaxation. In our model, besides the first HH 
subband, we also consider the contribution from the first 
LH subband and the second HH subband. Since only 
the hole states around the Fermi surface can contribute 
to the electron-hole Coulomb or exchange scattering, we 
choose d^^N\/ d^f^Nh as the criterion of the contribution 
from A hole subband. We further show the regime where 
the contribution from high hole subbands becomes signif- 
icant in Fig. [1] (the regime above the yellow curve), where 

c'Mh(^LH(i) + ^HH(2))/(^/ih^'i > O-l- It is noted that we 
only discuss the combined effect of the DP spin relaxation 
from the LH*^^^ and HH*^^) subbands in the following, as 
the effects on the electron-hole Coulomb scattering from 
these two subbands are analogous. Moreover, the ma- 
trix elements in Eq. ([5]) relevant to the HH^^^ subband 
are one order of magnitude smaller than those relevant 
to the LH*^^) subband for the relevant range of center-of- 
mass momentum K in the following cases. Therefore, we 
only discuss the effect on the BAP spin relaxation from 
the LH(^' subband. 

We first show how the multi-hole-subband effect in- 
fluences the temperature dependence of the spin relax- 
ation. The SRTs due to the DP and BAP mechanisms 
as well as the ratio tbap /tdp are plotted in Fig. [3] as 
function of temperature for a typical case with iV^ = 0, 
Nh — 5 X 10^^ cm~^ and A^ox = 10^^ cm~^. It is seen 
that after considering the contribution from high hole 
subbands, tbap decreases but tdp increases, and hence 
the importance of the BAP mechanism is enhanced. The 
underlying physics is as follows. The states in high hole 
subbands provide additional scattering channel, and the 
electron-hole Coulomb and exchange scatterings are both 
enhanced. The former suppresses the DP mechanism in 
the strong scattering limit, and the latter leads to an 
enhancement of the BAP mechanism. Both make the 
BAP mechanism become more important compared with 
the DP mechanism. It is also seen that the multi-hole- 
subband effect becomes more pronounced at high tem- 
perature. This is because the occupation of the high hole 
subbands becomes larger when temperature increases. 

From Fig. [21 one also finds that the multi-hole-subband 
effect does not significantly affect the trend of the tem- 
perature dependence of the SRT. The main change after 
the inclusion of the high hole subbands is that the tem- 
perature at which tbap /tdp reaches minimum becomes 
closer to the hole Fermi temperature. The underlying 
physics is as follows. In the degenerate regime (T < Tp), 
it is seen that compared with those in the single-hole- 
subband model, tdp (tbap) in the multi-hole-subband 
model increases (decreases) faster with increasing tem- 
perature, both originate from the increase in the occu- 
pation of the high hole subbands and hence the increase 
of the electron-hole Coulomb and exchange scatterings. 
This leads to a faster decrease of tbap /tdp with increas- 
ing temperature when T < Tp .Mi SS Nevertheless, in the 
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FIG. 3; (Color online) SRTs due to the DP and BAP mecha- 
nisms, the total SRT together with the ratio tbap/tdp vs. 
temperature T with A''i = 0, A^h = 5 x 10^^ cm~^ and 
Acx = 10^^ cm~^. The hole and electron Fermi temperatures 
are Tp = 124 K and Tp = 41 K, respectively. The sold curves 
are the results calculated with the lowest two subbands of HH 
and the lowest subband of LH. The dash curves are those from 
only the lowest subband of HH. Note the scale of tbap/tdp 
is on the right-hand side of the frame. 



nondegenerate regime (T > Tp), it is seen that tdp 
(tbap) in the multi-hole-subband model decreases faster 
(slower) than that in the single-hole-subband model. The 
accelerating in the decrease of tdp can be understood as 
follows. In the nondegenerate regime, the electron-hole 
Coulomb scattering decreases with temperature. With 
the contribution from high hole subbands, the electron- 
hole Coulomb scattering becomes stronger and thus the 
decrease rate also becomes larger. Therefore, tdf de- 
creases faster in the multi-hole-subband calculation i^^i^° 
The slowdown in the decrease of tbap originates from the 
anomalous K dependence of the matrix element M~i i 

2 ' 2 

in Eq. (O, which is relevant to the LH^^^ subband. As 
discussed above, in the nondegenerate regime, the tem- 
perature dependence of tbap is mainly from the matrix 
elements. It is also noted that the magnitude of M~i i 

2 ' 2 

decreases with K, whereas the magnitudes of the other 
matrix elements increase with K. With the increase of 
temperature, more holes and electrons are distributed 
at larger momentums, the contribution from M~i i de- 

2 ' 2 

creases, while the contributions from the other matrix el- 
ements increases. These two trends counteract with each 
other and make tbap decrease with increasing T very 
slowly at high temperature. Consequently, when T > Tp , 
the ratio tbap /tdp shows a steeper increase with the ris- 
ing temperature in the multi-hole-subband calculation. 
Therefore, both trends when the temperature is below 
and above Tp make the minimum of tbap /''dp appear at 
the temperature closer to Tp in the multi-hole-subband 
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FIG. 4: (Color online) SRTs due to the DP and BAP mecha- 
nisms, the total SRT together with the ratio tbap /tdp vs. 
hole density Nu with impurity densities Ni = Q (o) and 
Ni = Nh (■) at (a) T = 10 K and (b) 100 K (note that the 
scale of tbap/tdp is on the right-hand side of the frame). The 
excitation density Acx ~ 10® cm~'^. The sold curves are the 
results calculated with the lowest two subbands of HH and the 
lowest subband of LH. The dashed curves are those calculated 
with only the lowest subband of HH. The two purple dotted 
vertical lines indicate the hole densities satisfying i5p = ksT 
and Ep — AiJ^jjCi) = ^bT, respectively [note that the sec- 
ond dotted line in (b) is very close to the right frame]. Here 
A£'ljj(i) represents the energy splitting between the LH'^' 
and HH'^' subbands. 



calculation. 

We also investigate the multi-hole-subband effect on 
the hole-density dependence of the spin relaxation. In 
Fig. m the SRTs due to various mechanisms, the total 
SRT together with the ratio tbap /tdp are plotted as 
function of hole density. It is interesting to see from 
Fig. m^a) that at low temperature, the spin relaxation 
has a very intriguing hole-density dependence. In the 
impurity-free case, the hole-density dependence of the 
total SRT shows a double-peak structure: i.e., it first in- 
creases slightly, then decreases, again increases rapidly 
and finally decreases with increasing hole density. In 
the high impurity density case with Ni = the total 
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SRT shows first a peak and then a valley as a function 
of hole density. We first discuss the impurity-free case 
where the BAP mechanism is negligible and the double- 
peak structure is solely from the DP mechanism. The 
first peak can be understood as follows. The electron- 
hole Coulomb scattering increases with Nh in the non- 
degenerate regime (E'p < k^T) from the increase of the 
hole distribution, but decreases with Nh in the degener- 
ate regime (E'p > UbT) due to the increase of the Pauli 
blocking of holes i^^i^° Hence tqp first increases and then 
decreases with with the peak appearing around the 
hole density satisfying — k-oT. This behavior is sim- 
ilar to the peak predicted in n-type samples;^ where 
the peak originates from the electron-electron Coulomb 
scattering and hence appears around the electron den- 
sity satisfying = k-gT. It is also seen that the sec- 
ond peak only appears in the multi-hole-subband cal- 
culation, but becomes absent in the singlc-hole-subband 
calculation (the green dashed curve with circles), which 
indicates that this peak comes from the contribution of 
the electron-hole Coulomb scattering from high hole sub- 
bands. In fact, the scenario is similar to the first one. 
When Nh > 6 x 10"'^^ cm~^, the contribution from the 
LH^^) subband becomes important4^ Since the holes in 
the LH^^^ subband are still in the nondegenerate regime, 
the electron-hole Coulomb scattering increases with in- 
creasing hole density. Thus tdp increases with Nh ■ When 
TVh > 9 X 10" cm-2, i.e., E^ - Ai^LHd) > with 
A£^^Tj(i) representing the energy splitting between the 
LH^i) and HH^i) subbands, the holes in the LH^^) sub- 
band are also in the degenerate regime, and thus the 
effect of the Pauli blocking becomes significant. Conse- 
quently Tdp decreases with Nh- The second peak appears 
aroimd the hole density satisfying E^ — A£'ljj(i) = k-^T. 

Then we turn to the case of high impurity density with 
the impurity density being identical to the hole density. 
The scenario of the peak is as follows. The SRT due 
to the DP mechanism increases monotonically with hole 
density, since the electron-impurity scattering increases 
with Nh{= Ni). Moreover, the SRT due to the BAP 
mechanism decreases with Nh due to the increase of the 
hole distribution of the HH'^^^ and LH^^^ subbands, sim- 
ilar to the electron-hole Coulomb scattering. As a re- 
sult, the peak appears around the hole density where 
the BAP mechanism surpasses the DP mechanism.^*' It 
is also seen that Ttot increases with hole density when 
Nh > 9x 10^^ cm~^. The underlying physics is as follows. 
Similar to the previous discussion of the temperature de- 
pendence, with the increase of hole density, the decrease 
of the matrix element M~i i counteracts the increase 

2 ' 2 

of the other matrix elements. Thus the dependence of 
the hole density from the matrix elements is weak. Con- 
sequently, when the holes in the 1111*^^^ and LH^^^ sub- 
bands are both in the degenerate regime, the effect of 
the increase of the Pauli blocking is dominant, and tbap 
increases with Nh- The valley appears around the hole 
density satisfying Ep — AEj^^ii^) = ksT. It is also noted 
that the increase of tbap only appears in the multi-hole- 



subband calculation. In the single-hole-subband calcula- 
tion, Tbap does not increase with Nh but remains almost 
a constant for high hole density (the blue dashed curve) , 
since the increase of the Pauli blocking is counteracted 
by the increase of the matrix elements relevant to the 
HH(i) subband. 

The hole-density dependence of the spin relaxation at 
high temperature is also investigated [Fig. Sfb)] . Differ- 
ing from the behavior at low temperature, it is seen that 
there is only one peak in both the impurity-free and high 
impurity density cases. We further show that the peaks 
in both cases come from the competition of the DP and 
BAP mechanisms, which is similar to the peak in the 
case with high impurity density and low temperature. 
The absence of the peak from the electron-hole Coulomb 
scattering is due to the multi-hole-subband effect. As 
discussed above, when the hole density is high enough so 
that the holes in the lowest subband are in the degener- 
ate regime, the contribution of the electron-hole Coulomb 
scattering from the lowest hole subband decreases with 
increasing hole density due to the increase of the Pauli 
blocking. However, at high temperature the contribution 
from high hole subbands is also important in this hole 
density regime. It is further noted that the holes in the 
high subbands are still in the nondegenerate regime, thus 
the contribution from the high hole subbands increases 
rapidly with Nh and surpasses the effect from the lowest 
hole subband. Consequently, tdp continues to increase 
with Nh and the Coulomb peak disappears. 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have performed a comprehensive in- 
vestigation of electron spin relaxation in p-type GaAs 
quantum wells from a fully microscopic KSBE approach. 
All relevant scatterings, such as, the electron-impurity, 
electron-phonon, electron-electron Coulomb, electron- 
hole Coulomb, and electron-hole exchange (the BAP 
mechanism) scatterings are explicitly included. 

We present a phase-diagram-like picture showing the 
parameter regime where the DP or BAP mechanism is 
more important. In the case with no impurity and high 
excitation density, our results are consistent with those in 
Ref. [T^: i.e., the BAP mechanism is unimportant at low 
temperature, which is in stark contrast with the common 
belief in the literaturciiii'^'^i^ However, since we extend 
the scope of our investigation to higher hole density by in- 
cluding more hole subbands in the model, it is discovered 
that the BAP mechanism can surpass the DP mechanism 
in the regime with high temperature and sufficiently high 
hole density. In the cases with low excitation density 
and/or high impurity density, the regime where the BAP 
mechanism surpasses the DP mechanism becomes larger. 
We also show that the temperature regime where the 
BAP mechanism is more efficient than the DP mechanism 
is very different in the impurity-free and high impurity 
density cases. In the impurity-free case, this regime is 
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around the hole Fermi temperature for high hole den- 
sity, regardless of excitation density. However, in the 
high impurity density case with the identical hole and 
impurity densities, this regime is roughly from the elec- 
tron Fermi temperature to the hole Fermi temperature. 
This is because the dominant scatterings in these two 
cases are the electron-hole Coulomb scattering and the 
electron-impurity scattering, respectively. The different 
dominant scatterings lead to the different temperature 
dependences of tdp, and hence the different behaviors 
of the ratio tbap/'^dp- In particular, in the case with 
high impurity and very low electron excitation densities, 
the electron (hole) Fermi temperature is much lower than 
(close to) the lowest (highest) temperature of our inves- 
tigation. Consequently, the BAP mechanism can domi- 
nate the spin relaxation in the whole temperature regime. 
Moreover, we predict that for the impurity-free case, in 
the regime where the DP mechanism dominates the spin 
relaxation, e.g., the cases with high excitation or low hole 
density, the total SRT presents a peak around the hole 
Fermi temperature, which is from the nonmonotonic tem- 
perature dependence of the electron-hole Coulomb scat- 
tering. 

The multi-hole-subband effect on the spin relaxation is 
also revealed. It is shown that the multi-holc-subband ef- 
fect enhances the relative importance of the BAP mecha- 
nism significantly for high temperature and/or high hole 
density. We also predict that at low temperature the 
spin relaxation has a very intriguing hole-density depen- 
dence thanks to the contribution from high hole sub- 
bands. In the impurity-free case, the total SRT shows 
a double-peak structure. Both peaks originate from the 
nonmonotonic hole density dependence of the electron- 
hole Coulomb scattering. The only difference is that the 
first (second) peak comes from the contribution from the 
HH(i) (LH(i)) subband, and hence appears around the 
hole density satisfying = ksT {E^-AE^^n) = k^T). 
In the high impurity density case with identical impu- 
rity and hole densities, there are first a peak and then a 
valley. The peak is formed as the DP and BAP mecha- 
nisms compete with each other: the SRT due to the DP 
(BAP) mechanism increases (decreases) with T as the 
DP (BAP) mechanism dominates at low (high) tempera- 
ture. Moreover, since the decrease of the matrix element 
M~i I counteracts the increase of the other matrix el- 

2 ' 2 

ements of the BAP scattering, the hole-density depen- 
dence from the matrix elements is weak. Consequently, 
when the holes in the HH*^^) and LH^^^ subbands are both 
in the degenerate regime, the effect of the increase of 
the Pauli blocking is dominant, and tbap increases with 
Nh- Therefore the valley is formed. However, at high 
temperature, we show that the peak from the electron- 
hole Coulomb scattering disappears and only the peak 
from the competition of the BAP and DP mechanisms 
remains. 

Finally, we discuss the relevance to experiments. Since 
electrons are minority carriers in p-type semiconduc- 
tors, they have finite lifetime limited by the electron- 



hole recombination. From experiments, the electron life- 
time is found to be on the order of several hundreds 
of picoseconds which may be shorter than the spin 
relaxation time in some parameter regime such as the 
high impurity density case. However, it has been demon- 
strated that spin relaxation time much longer than the 
lifetime can be measured via the Hanle or other time- 
resolved measurements fii^ Hence the calculated long 
spin relaxation times can be observed in reality. In short, 
our calculation based on realistic parameters provides 
useful information in a wide range of experimental con- 
ditions, which would benefit the understanding on spin 
relaxation. 
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FIG. 5: (Color online) SRTs r from the experimental data 
in Ref. and from the DP and BAP mechanisms calculated 
via the KSBE approach in a p-type GaAs quantum well. 



APPENDIX: COMPARISON WITH 
EXPERIMENT 

We compare the calculation via the KSBE approach 
with the experimental results in Ref. 50. In Fig. [5l we 
plot the temperature dependence of the SRTs from our 
computation and from the experiment in p-type GaAs 
quantum well. Here a = 3 nm, Nh = 3 x 10^^ cm^^ and 
A'ex = 2.6 X 10^ cm^^ as indicated in the experiment.— 
Ni — 0.3Nh is obtained by fitting the mobility /i ~ 
4000 cm^V-^s"^ given in Ref. [H We take the Dressel- 
haus spin-orbit coupling parameter as 7d = 18.3 eV-A^. 
As the range of 7d in bulk GaAs calculated and measured 
via various methods is from 6.4 to 25.5 eV-A^j^ our value 
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of 7d is reasonable. It is seen that our calculation is in 
good agreement with the experimental data. The devi- 
ation at T = 5 K is likely to be from the fact that the 
electron temperature is higher than the lattice tem- 
perature T in experiment due to photo-excitation. This 
effect becomes significant at low temperature.— It is also 



noted that the SRT due to the BAP mechanism is over 
one order of magnitude larger than that from the DP 
one, i.e., the spin relaxation is totally governed by the 
DP mechanism in this case. This is consistent with our 
conclusion that the BAP mechanism is unimportant at 
low impurity and hole densities. 



* Author to whom correspondence should be addressed; 
Electronic address: mwwu@ustc.edu.cn. 

^ F. Meier and B. P. Zakharchenya, Optical Orientation 
(North-Holland, Amsterdam, 1984). 

^ S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 
Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka- 
nova, and D. M. Treger, Science 294, 1488 (2001). 

^ Semiconductor Spintronics and Quantum Computation, 
edited by D. D. Awschalom, D. Loss, and N. Samarth 
(Springer- Verlag, Berlin, 2002); I. Zutic, J. Fabian, and 
S. Das Sarma, Rev. Mod. Phys. 76, 323 (2004); J. Fabian, 
A. Matos-Abiague, C. Ertler, P. Stano, and I. Zutic, Acta 
Phys. Slov. 57, 565 (2007); Spin Physics in Semiconduc- 
tors, edited by M. I. D'yakonov (Springer, Berlin, 2008); 
and references therein. 

* G. L. Bir, A. G. Aronv, and G. E. Pikus, Zh. Eksp. Teor. 
Fiz. 69, 1382 (1975) [Sov. Phys. JETP 42, 705 (1976)]. 

^ M. I. D'yakonov and V. I. Perel', Zh. Eksp. Teor. Fiz. 60, 

1954 (1971) [Sov. Phys. JETP 33, 1053 (1971)]; Fiz. Tverd. 

Tela (Leningrad) 13, 3581 (1971) [Sov. Phys. Solid State 

13, 3023 (1972)]. 
^ M. W. Wu and C. Z. Ning, Eur. Phys. J. B 18, 373 (2000); 

M. W. Wu and H. Metiu, Phys. Rev. B 61, 2945 (2000); 

M. W. Wu, J. Phys. Soc. Jpn. 70, 2195 (2001). 
'' A. G. Aronov, G. E. Pikus, and A. N. Titkov, Zh. Eksp. 

Teor. Fiz. 84, 1170 (1983) [Sov. Phys. JETP 57, 680 

(1983)]. 

® G. Fishman and G. Lampel, Phys. Rev. B 16, 820 (1977). 
^ K. Zerrouati, F. Fabre, G. Bacquet, J. Bandet, J. Frandon, 
G. Lampel, and D. Paget, Phys. Rev. B 37, 1334 (1988). 
^° P. H. Song and K. W. Kim, Phys. Rev. B 66, 035207 

(2002) . 

" M. Z. Maialle, Phys. Rev. B 54, 1967 (1996). 

^2 J. Zhou and M. W. Wu, Phys. Rev. B 77, 075318 (2008). 
For brief review, see M. W. Wu, M. Q. Weng, and J. L. 
Cheng, in Physics, Chemistry and Application of Nanos- 
tructures: Reviews and Short Notes to Nanomeeting 2007, 
eds. V. E. Borisenko, V. S. Gurin, and S. V. Gaponenko 
(World Scientific, Singapore, 2007) p. 14, and references 
therein. 

" M. Q. Weng and M. W. Wu, Phys. Rev. B 68, 075312 

(2003) ; ibid. 70, 195318 (2004). 

M. Q. Weng, M. W. Wu, and L. Jiang, Phys. Rev. B 69, 
245320 (2004). 

L. Jiang and M. W. Wu, Phys. Rev. B 72, 033311 (2005). 
J. Zhou, J. L. Cheng, and M. W. Wu, Phys. Rev. B 75, 
045305 (2007). 

^® P. Zhang, J. Zhou, and M. W. Wu, Phys. Rev. B 77, 
235323 (2008). 

^3 J. H. Jiang and M. W. Wu, Phys. Rev. B 79, 125206 
(2009). 

2° J. H. Jiang, Y. Zhou, T. Korn, C. Schiiller, and M. W. Wu, 
Phys. Rev. B 79, 155201 (2009). 



X. Z. Ruan, H. H. Luo, Y. Ji, Z. Y. Xu, and V. Umansky, 

Phys. Rev. B 77, 193307 (2008). 
^2 L. H. Teng, P. Zhang, T. S. Lai, and M. W. Wu, Europhys. 

Lett. 84, 27006 (2008). 
2^ D. Stich, J. Zhou, T. Korn, R. Schulz, D. Schuh, W. 

Wegscheider, M. W. Wu, and C. Schiiller, Phys. Rev. Lett. 

98, 176401 (2007); Phys. Rev. B 76, 205301 (2007). 

24 T. Korn, D. Stich, R. Schulz, D. Schuh, W. Wegscheider, 
and C. Schiiller, Adv. Solid State Phys. 48, 143 (2009). 

25 D. Stich, J. H. Jiang, T. Korn, R. Schulz, D. Schuh, W. 
Wegscheider, M. W. Wu, and C. Schiiller, Phys. Rev. B 
76, 073309 (2007). 

A. W. HoUeitner, V. Sih, R. C. Myers, A. C. Gossard, and 
D. D. Awschalom, New J. Phys. 9, 342 (2007). 

F. Zhang, H. Z. Zheng, Y. Ji, J. Liu, and G. R. Li, Euro- 
phys. Lett. 83, 47006 (2008); ibid. 83, 47007 (2008). 

2* M. KrauB, R. Bratschitsch, Z. Chen, S. T. Cundiflt, and H. 

C. Schneider, larXiv:0902.0270 
2^ C. Yang, X. Cui, S.-Q. Shen, Z. Xu, and W. Ge, Phys. 

Rev. B 80, 035313 (2009). 
^° K. Shen, Chin. Phys. Lett. 26, 067201 (2009). 
^1 E. T. Yu, J. O. McCaldin, and T. C. McGill, Solid State 

Phys. 46, 1 (1992). 
^2 H. Haug and A. -P. Jauho, Quantum kinetics in Transport 

and Optics of Semiconductors (Springer, Berlin, 1998). 

G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 
(1984); JETP Lett. 39, 78 (1984). 
^5 W. H. Lau and M. E. Flatte, Phys. Rev. B 72, 161311(R) 
(2005). 

^'^ A. N. Chantis, M. van Schilfgaarde, and T. Kotani, Phys. 
Rev. Lett. 96, 086405 (2006). 

G. D. Mahan, Many Particle Physics (Kluwer Academic, 
New York, 2000). 

M. Z. Maialle, E. A. de Andrada e Silva, and L. J. Sham, 
Phys. Rev. B 47, 15776 (1993). 

The contribution from the short-range term is negligibleii^ 
*° Our calculation is not sensitive to the initial spin polariza- 
tion when it is small. 
4^ Semiconductors, Landolt-Bornstein, New Series, edited by 

O. Madelung (Springer, Berlin, 1987) Vol. 17a. 
''2 W. Ekardt, K. Losch, and D. Bimberg, Phys. Rev. B 20, 
3303 (1979). 

Most cases we discuss in this paper are in the strong scat- 
tering regime, where tdp increases with momentum scat- 
tering strength and decreases with inhomogeneous broad- 
ening. 

^'^ The hole density Nh consists of both holes from acceptors 
and thoes from optical excitation. 

45 It is seen that in the case with Aex = 10^^ cm"^^ the 
variation of the SRT due to the DP mechanism becomes 
very mild when T < 10 K. This is because this case belongs 
to the intermediate scattering regime, where tdp has a 



10 



weak dependence of the scattering strength. 

G. F. Giuhanni and G. Vignale, Quantum Theory of the 

Electron Liquid (Cambridge University Press, Cambridge, 

2005). 

M. M. Glazov and E. L. Ivchenko, Pis'ma Zh. Eksp. Teor. 
Fiz. 75, 476 (2002) [JETP Lett. 75, 403 (2002)]; Zh. Eksp. 
Teor. Fiz. 126, 1465 (2004) [JETP 99, 1279 (2004)]. 
F. X. Bronold, A. Saxena, and D. L. Smith, Phys. Rev. B 
70, 245210 (2004). 



It is also noted that the contribution from the HH^^' sub- 
band is always negligible in this case. 

M. D. Martin, L. Vina, M. Potcrnski, and K. H. Ploog, 
Phys. Stat. Sol. (b) 215, 229 (1999); L.Vina, J. Phys.: 
Condens. Matter 11, 5929 (1999). 

J. Wagner, H. Schneider, D. Richards, A. Fischer and K. 
Ploog, Phys. Rev. B 47, 4786 (1993). 



